On stability of a new model of wormhole 
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Abstract 



We investigate the stability of the wormhole of the Mors-Thorne-Ellis-Bronnikov 
type. In our model the matter inside it consists of a radial magnetic field and an 
ideal phantom- fluid. Properties of the matter are described in section H] HT1 and IIIII of 
this paper. We consider spherical perturbations only and find examples of the stable 
wormholes against these perturbations. 

I. INTRODUCTION 

The problem of the stability of the wormholes has been discussed in many papers (see [lj- 

id). 
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- |l3| (phantom scalar field with the negative kinetic 
[l, 2]) and in accordance with more 



The simplest model of a wormhole 
term) proved to be unstable (in spite of erroneous papers 

nnn 

recent (independent) correct research [3H7I. Il4j|. 

A more complex model of a wormhole (with the same metric 0, H| — radial magnetic 
field and the phantom dust with the negative energy density) is almost stable under all types 
of the spherical perturbation (except for the longitudinal radial motion of the dust on the 
inertia). The growth of unstable type is sufficiently slow: it increases proportionally to the 
time [8|, |9j. Therefore, the authors of this paper it has been suggested that the unstable 
type can be easily crushed by the introduction into the model additional parameters. In the 
paper [if]] it was found a stable thin-shell traversable wormholes model, with a thin shell of 
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phantom matter on its throat (by using the cut-and-paste procedure). 

All of these studies provided hope for the existence and finding fully stable solution (under 
all spherical perturbation types). 

In this paper we investigate the stability (under the spherical perturbations) the similar 
model of the Mors-Thorne-Ellis-Bronnikov type wormhole (MTEB), see also 15 1. 

The matter of our model of the wormhole (WH) consists of the radial monopole magnetic 
field (with topological charge q) and the ideal phantom fluid. In the stationary case (respect 
to which will be a study on stability) the energy density e of this fluid is negative and equal 
to twice of the absolute value of the energy density of the magnetic field. The pressure p of 
the fluid in a stationary case is equal to zero (phantom dust). In the case of deviation of the 
energy density e from its stationary value — the pressure p is proportional to this deviation 
/ oc e — Bp- As already mentioned, without pressure, this model is linearly unstable over 
time |g, 9|. 



II. EQUATIONS OF THE MODEL 

It is convenient to choose the metric tensor in the spherically symmetric case like the 
following 1 : 

ds 2 = e u dt 2 -e x dx 2 -e a {dd 2 + sin 2 dip 2 ). (1) 

Here e a = r 2 , where 4irr 2 - area of a sphere around the center of the system, r, v and A are 
functions x and t. 

The Einstein equations corresponding to the metric in the comoving reference sys- 



Units are chosen as: c = 1 and G = 1 (the speed of light and gravitational constant). 
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tern 17j can be written as 2 : 



lire + q 2 e~ 2a 



a" + -(a') 
4 
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Inp - q 2 e~ 2a = -e~ v 
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lirp + q 2 e~ 2a = e~ x 



-v" + -{V 
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2a' + da' — Ac/ — v 'a = . 



(4) 
(5) 



The stationary wormhole MTEB metric tensor is defined by 



ds 2 = dt 2 -dx 2 -(q 2 + x 2 ) ■ (d9 2 + sin 2 6 dip 2 ). 



(6) 



This WH has the energy- momentum tensor which corresponds to a mixture of monopole 
electric (or magnetic) field and the phantom dust: 



87rT™ = diag 



+q 2 



+q 2 



(g 2 + x 2 ) 2 ' (g 2 + x 2 ) 2 ' (g 2 + x 2 ) 2 ' (q 2 + x 2 'f 



+ diag 



-2q 2 



(7) 



The first term on the right side corresponds to the energy-momentum tensor of the electric 
(or magnetic) field with a charge q; second - dust matter with negative energy density. 

III. LINEARIZATION OF THE EQUATIONS 



Let us consider small spherical perturbations of matter and metrics of the MTEB worm- 
hole. 

We introduce the notations: 
-2g 2 



lire 



+ f(x,t), e a = (q 2 + x 2 )e^\ £ 2 ^(g 2 + x 2 ), Srcp = hf. (8) 



(q 2 + x 2 ) 2 

Here h(x) — arbitrary function, its physical meaning — a square of sound speed v 2 (in this 
fluid). 

We introduce the dimensionless coordinates: namely we put q = 1. 



2 The derivation of these equations can be found, for example in [ItJi (task 5 in §100), prime denotes the 
derivative with respect to x, and dot - with respect to t. 
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We write the equations (EH) in the linear approximation for small perturbations: 
v, \ V, f- 



A - rj - 3xrj' + xX' „ 2rj + \ 
f = 72 V + (9) 



, „ ri(l — x 2 ) Xx 2 xiv' + ??') , . 

hf + v+ 74 72 = °> ( 10 ) 



v" + rj" — A — 77 , 2ry — A x(rj' - X' /2 + u' /2) 

I 2 

f)' + x{r)-X)/£ 2 = 0. (12) 



— o + t4 + t2 ' v 11 ' 1 



From ( II 2p we have: 



e 2 

v - X = Fi(x) - — rf (13) 



If the pressure p is isotropic, then from equations (J2]|5]) we can obtain two useful relations 
which derived directly from the formula T^ k = (consequence of Bianchi identities): 

X + 2a = -2e/(p + e) , v' = -2p'/(p + e) . (14) 

Or in the linear approximation: 

A + 2f, = ef, v' = i\hf)'. (15) 

From (II 5p we have: 

e 4 / = X + 2r ] + F 2 (x) (16) 
Expressing A from (fTBl . and taking into account (ITB"]) from (Q we have: 

F 2 (x) = (17) 

The arbitrary functions Fi(x) and F 2 (x) are determined by the choice of initial conditions 
(at t = 0) to small perturbations A, / and v (for given initial condition at r)(x, t)). Therefore, 
renaming at t = in f[T3"j) and f[T6"j) for small perturbations: A — > A — Fi and /—>•/ + F 2 /£ t A 
then we can put the function Fi and F2 equal to zero. 
From ([9]), ([TOD and (II3D we have: 

/ = 3r ? /e 4 + V/K 2 )= [Z 3 v]'/te% (18) 
^ + /i/-xe 2 (/i/)' + r / /e 4 = 0, (19) 
A = 1 [ft]' • (20) 

37 



5 



From the equations ( |T8lfT9l) we obtain: 



rj — hrj" + 7/ 



where 



2h 

x£ 2 



-ti 



+ rflJ(x) = 



(21) 



U(x) = 



12hx 2 + 3h - 3tix£ 2 + 1 



(22) 



Equation (|21|) at h > is a differential equation of hyperbolic type. At infinity (x — > ±00) 
the potential U(x) tends to zero and equation (121]) becomes an equation for sound waves 
with the speed of sound v s = \fh. 



IV. STABILITY INVESTIGATION 

We transform the equation fl2TT) to the canonical form. To do this we change variables 
x — > z. Then d x — > f3d z , where j3 = 1^. Equation ([21]) can be rewritten as: 



77 - 



2/2 

/i/9 2 r?, 22 -/3r/, z ( — - - ti 



xC, 



+ rjU = (23) 



Now choose a function h(x) and /3(x), so that in the square brackets expressions (1231) was 

only value of 77, zz , ie /i/3 2 = 1 and h(3' — /3[2/i/(x£ 2 ) — ti] = 0. 

From these conditions we obtain h = hox A /^ A , where ho = const > 0. 

Then the hyperbolic equation (|2TT) can be rewritten in the canonical form: 

V-V,zz+vU(z)=0 (24) 

where 



12h x 6 - 9h x A + (1 + x 2 ^ 2 

IT+x~ 2 



U \ Z ) = H , ( 25 ) 



Here we took into account that now the variable x depends from the variable z by the 
following equation: 

x 2 - y/h^zx -1 = (26) 

and 



_,_ _ y/h~o~z ± Vhoz 2 + 4 



= ^ ^ (27) 
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Here the sign "+" corresponds to the x > 0, and the sign "— " - - respectively x < 0. We 
first investigate the region x > 0. 

Asymptotics of z — > — oo corresponds to x — > (y/h ~\z\)~ 1 — > 0, and asymptotics z — > +oo 
corresponds to x — > {^/Tlqz) — > +oo. 

The solution of equation (|24|) can be obtained by separation of variables: 

oo 

V (z,t) = J2Tn(m n (z), (28) 

n=0 

tT = ~© U(z) = -w n . (29) 

From ff29l we obtain: T n = exp(iw n t). Here, the quantity w n = const has a physical sense 
to the harmonic oscillation frequency with the number n, for the small perturbation rj. 
From ( 1291) we obtain for each harmonic: 

* n>zz + w 2 J? n - U(z)V n = (30) 

The expression flBDl) is a stationary Schrodinger equation for the whole numerical axis 
(— oo, +oo) with the potential U(z). 



As is well known (see. 18j . §18) energy levels E n = (w n h) 2 /2 in the Schrodinger operator 
spectrum are always positive, if the effective potential U(z) is regular and non-negative (for 
all z) and U(z) — > at z — > oo. From ( 1241) follows that U(z) satisfies these conditions at 
< h < h 00 2.8 — see figure. However, the physical constraint on the parameter ho is 
the maximum speed of sound, which can not exceed the speed of light. Since the square of 
the speed of sound is || = h < c 2 , then there must be h max = h < 1. 

The result is w 2 > 0, ie oscillation frequency must be real value. 
The proof for region x < is analogous. 

Thus we demonstrated the existence of the matter model in which the function r)(x,t) is 
a nonincreasing. 

Since each term of ( 1281) nonincreasing for time evolution, there must be a non-increasing, 
and the derivative by x from each member of this series (as derivative by x does not affect on 
the time component T n , which responsible for the time dependence). Therefore, we can make 
a statement that the quantity r)'(x,t) well as nonincreasing (at < h <l). Therefore, the 
functions u(x,t), f(x,t) and \(x,t) are nonincreasing (according to expressions ( fT5l [T8H201 ). 
Thus, the complete solution is stable. 
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(a) Side view (b) Top view 

FIG. 1: View of the surface for the potential U(x,ho) for h[x) = Hqx 4 /^ 4 . This shows that region 
U > corresponds to the range < ho < /loo ~ 2.8 for any values of x. 

The case h = should be considered individually: Eq. (123]) at h — becomes an equation 
of parabolic type and the variables (x, t) are not separated, as decomposition in the form 
fl28|) becomes inapplicable. The solution of equation ( 12 ip with h = is easily found in the 
form: r] h=0 (x,t) = r] exp (it/£ 2 ). This solution is also stable. However, the derivative of 
this solution r]' h=0 = —2itxr] h=0 /^ 4 at x ^ is not asymptotically stable function, and has a 
linear instability with respect to time. Consequently, the functions Xh=o(x, t) and fh=o(x,t) 
at x 7^ are also linearly unstable (see. Q, 1^]). 



V. CONCLUSIONS 



We have demonstrated that in general relativity, it is possible to construct a model of static 
and traversable wormhole, which will be stable with respect to small spherical perturbations. 

In this paper we do not perform investigation on the stability by non-spherical perturba- 
tion types, but it is known (see. 19j), that the non-spherical perturbation types, seems to be 
more stable than the spherical, because they have a centrifugal (and other higher-multipole) 
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barriers in the effective potential for perturbations. 
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